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1. Introduction 

Pions are the lightest hadrons in nature. Their special status as the lightest members of an 
octet of pseudoscalar (pseudo-)Goldstone bosons associated with the breaking of the SU (3) a part 
of chiral symmetry makes them fascinating objects to study. In effective theories, pions appear as 
elementary degrees of freedom. Naturally, these theories have nothing to say about the internal 
structure of the pion. Within QCD, however, pions are formed as bound states of a quark and 
an anti-quark interacting strongly with each other. Powerful constraints such as the axial Ward- 
Takahashi identity guarantee the massless Goldstone boson nature of the pions in the chiral limit, 
although the constituents remain massive. This fascinating interplay between dynamical mass gen- 
eration on the quark level and the formation of massless bound states is particularly transparent 
in functional approaches to QCD using Dyson-Schwinger, Bethe-Salpeter and Faddeev equations, 
see e.g. [1-4] for reviews. 

An important probe for the internal structure of hadrons is their response to electromagnetic 
fields. Real and virtual photons serve to extract quantities such as form factors, polarizabilities and 
distribution functions which in turn offer interesting insights into global as well as spatially resolved 
electromagnetic properties of these hadrons. Virtual pions (the 'pion cloud' [5]) play an important 
role in this respect. At small momentum transfer, the incoming photons couple dominantly to the 
virtual pion cloud of the bound state instead of the quark and antiquark constituents. An impressive 
demonstration of this effect has been given in [6], where the electromagnetic form factors of the 
nucleon have been determined in a three-body Faddeev approach. A similar picture has been found 
for the axial and pseudoscalar form factors [7]. The inclusion of pion cloud effects into form factor 
calculations within the functional approach is challenging, both conceptually and numerically. First 
steps in this direction have been made in [8, 9], where pion cloud effects in the masses of light 
mesons have been studied. In this proceedings contribution we report on our efforts to determine 
meson form factors in this approach. 

Another interesting place where virtual pions may play a role are tetraquark bound states. The 
story of the tetraquark is tightly connected to the issue of the light scalars and especially to the 
case of the a meson. Back in the 70ies, Jaffe introduced a simple quark-bag model to describe 
these states [10]. He showed that tetraquarks allow for a natural explanation for the inverted mass 
spectra and the huge width found in the lightest + meson nonet. Unfortunately this huge width 
and the multitude of other particles with + quantum numbers in the low mass region made a clear 
experimental signal difficult to obtain. Even the existence as a bound state, described as a pole in 
the S-matrix, was debated for decades. Only recently new experimental data coming from BES 
[11] and KLOE [12] combined with the application of a modified Roy-equation [13, 14] technique 
asserted the existence of the a meson. 

The mass of the a// (500) was found to be m a = 450 + /280MeV [13, 14] and the particle 
has been promoted to a proper member of the PDG again [15]. While the existence of the /o(500) 
seems beyond doubt, its interpretation in terms of its internal structure is still a matter of debate. 
Several approaches, including effective theory and large 7Y C studies [16-21] as well as more recent 
lattice calculations [22, 23], indicate a strong contribution from a qqqq-like. operator. 

In this contribution we report on our study of tetraquarks within an approach using a covari- 
ant four-body equation, which we approximate by a diquark and meson constituent picture. These 
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Figure 1: The Dyson-Schwinger equation for the fully dressed quark propagator. 



are determined self-consistenly from the the underlying quark and gluon substructure employing a 
model for the quark-gluon interaction known to reproduce hadron properties on a phenomenologi- 
cal level [24]. 

2. Meson form factors 

2.1 Dyson-Schwinger and Bethe-Salpeter equations beyond rainbow-ladder 

The Dyson-Schwinger equation (DSE) for the quark propagator is given diagrammatically in 
Fig. 1. On the left-hand side we find the fully dressed inverse quark propagator S^ 1 (p) = ipA(p ) + 
B(p 2 ). The first term on the right-hand side denotes its bare counterpart Sq 1 (p) = if + m, and the 
self-energy includes again the fully dressed quark propagator as well as the gluon propagator and 
a bare and a dressed quark-gluon vertex. 

For phenomenological purposes, a commonly used approximation scheme for the coupled 
system of the quark DSE and Bethe-Salpeter equation (BSE) is the rainbow-ladder approach. In 
the quark DSE this approximation amounts to the simple choice of T(p, q)fj. = Y^i r(&) for the quark- 
gluon vertex, where p and q are the ingoing and outgoing quark momenta and k is the momentum 
of the gluon. This approximation ignores further tensor structures of the full vertex as well as the 
dependence of the dressing function on the quark momenta. The remaining dressing T(k) can then 
be combined with the dressing Z(k) of the Landau-gauge gluon propagator, 



into a single effective coupling for the quark-gluon interaction. The most important property of this 
approximation is that it readily allows for the construction of a corresponding interaction kernel in 
the Bethe-Salpeter equation of mesons, thus satisfying the axial-vector Ward-Takahashi identity 
(axWTI). This is mandatory to obtain the pion in the chiral limit as both, a Goldstone boson and a 
bound state of a massive quark-antiquark pair. 

Although there is much phenomenological success in describing mass spectroscopy of light 
hadrons, decay constants [25-27] and dynamical properties like meson form factors [28] with 
good agreement to experiment [29], the rainbow-ladder scheme is not sufficient when it comes 
to the inclusion of unquenching effects like pion-cloud contributions. In the tower of coupled 
Dyson-Schwinger equations, such contributions appear already in the DSE for the quark-gluon 
vertex, see [9] for details. When these contributions are taken into account explicitly, the resulting 
quark-gluon interaction can be split into a pure Yang-Mills part (present already in the quenched 
approximation) and a separate part denoting the effects of the back-reaction of pions onto the quark 
propagator. The DSE is then given diagrammatically in Fig. 2 [9]. The corresponding equation is 




(2.1) 
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Figure 2: Dyson-Schwinger equation for quark propagator with explicit diagram for the pion back-reaction. 



given by 



s- 1 (p) = z 2 5 - 1 ( P ) + g 2 c F (z 2 ) 2 f ^ y,s{k) 7v (5, v - ^) 

d ' k l Z 2 y 5 S(k)rJ I ^-,k-p)+Z 2 y 5 S(k)rJ^-,p-k 



q^q v \ Z(q 2 )r YM (q 2 ) 



(2n 



D % {q 2 



(2.2) 



with the momentum routing q = p — k and the pion propagator D n (q 2 ) = 1 / (q 2 + M 2 ) . The product 
Z(k 2 )TyM{k 2 ) stands for the Yang-Mills part of the quark-gluon interaction, which we approxi- 
mate in the spirit of the rainbow-ladder approach. The Maris-Tandy model [27] that parametrizes 
Z(k 2 )TYM{k 2 ) is given by: 



471/71 

u ,6 



Z(k 2 )r YM (k z ) = I ^ Dk A exp{-k i /w z ) + 



2 2 27ry m [1 -exp(-<? 2 / 4m f 2 )] 



Mt + (1+^/A^ d )) 



(2.3) 



with m t = 0.5 GeV, z = e 2 - 1, y m = 12/(33 - 2%) and A QCD = 0.234GeV. In the pion part of the 
interaction, T n {p;P) is the full pion wave function which we approximate by its leading amplitude 
in the chiral limit [25]: 



ri{p;P) = zJy 5 



i B x(P 



(2.4) 



Here, B x (p 2 ) is the scalar dressing function of the quark propagator in the chiral limit. The re- 
placement of the leading physical pion amplitude by this chiral-limit approximation is correct on 
the few-percent level. Note that we use this approximation only for the internal pion which medi- 
ates the interaction. 

The Bethe-Salpeter equation for meson bound states is shown diagrammatically in Fig. 3. The 
corresponding equation reads: 



d A k 



Kl^(p,k;P)+K^ s (p,k;P)) \s(k + )r^(k;P)S(k_) 



(2.5) 



TV 






Figure 3: Bethe-Salpeter equation for quark anti -quark bound state with explicit diagram for the pion back- 
reaction. 
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with the kernels K™. s and Kf*°" s given by: 

<?> ,* ; P) = ^ 2 >r"V)(^ 



(s Wv \ 


\X a ' 




\X a ' 




k 2 ) 


Y 7,_ 


ts 


Y 7v _ 


ru 



q 2 

K?:° r n s (p,k;P) = l -[Yi] ru ^ P + k 2 ~ P ;p-k)[Z 2 ^Y 5 ] ts D 7l (q^ 

l rWl fp + k-P 
+ 4 



(2.6) 



rlU ( P + l~ P ' k -p) y 5 ] ts D n (q 2 ) (2.7) 

Here, T^{p;P) is the Bethe-Salpeter vertex function of a quark-antiquark bound state. The Latin 
indexes (t,u,r,s) of the kernels refer to their Dirac structure. It has been shown explicitly in [8] 
that this interaction kernel satisfies the axWTI. 

2.2 Coupling the meson bound state to an external field 

A systematic procedure to couple bound states to external gauge fields was given by [30] and 
applied to the diquark-quark model of baryons in Ref . [3 1 ] and to the three-body Faddeev equation 
in Refs. [6, 32]. In short, the evolution of the two-body quark system is given by the amputated 
version of the scattering matrix T. This function can be obtained by solving a Dyson equation: 

T = -iK-iKG T, (2.8) 

where Go is the disconnected product of two full quark propagators and —iK is the two-quark 
interaction kernel. When the two-quark system forms a bound state, the scattering matrix develops 
a pole at P 2 = —M 2 , and can be defined as: 

rr 

(2.9) 



P 2 +M 2 ' 

Substituting (2.9) in (2.8) and keeping only the singular term, we arrive at the Bethe-Salpeter 
equation for the two-quark bound state: 

r = -iKG r. (2.10) 
Then a systematic coupling to the external gauge field gives for T^: 

7 M = T{iK~ x K^K~ l +Gq)T . (2.11) 
The bound-state electromagnetic current J 11 can be expressed at the pole by: 

T»n 7 Tf 9 f ' (2.12) 

Pj+Mj Pf+Mj 

Substituting this in (2.1 1) and using (2.10), we get: 

j"=Tf(-iGoK^Go + G^)T i . (2.13) 
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So far, in rainbow-ladder approximation the first term Kg L = = because the gluon does not 

couple to a photon. Including the pion back-coupling, however, the gauged kernel does contribute, 
since the photon can couple to the exchanged pion. This fact generates two additional diagrams 
for the pion form factor, which are given in Fig. 4. Here is an ansatz for the pion-quark-quark- 
photon- vertex, which can be built along the corresponding one for the diquark-quark-quark-photon 
vertex derived in [31]. It reads: 

M " = qq A{v-q)-Q% iniP ~ q) ~ m ~ " 

+q - Z P -^Q Q) Q 2 {T{{P ~ q) + Q,1) - T{p - q)) - 

In comparison to rainbow-ladder, the calculation becomes more complicated. For example, the 
second diagram involves the pion form factor itself, so that will lead to the necessity to perform a 
self-consistent, iterative calculation which is also complicated by the two-loop integration. 

Up to now we have performed only a part of the calculation, namely the first diagram in Fig. 4. 
The resulting pion form factor with pion back-coupling is shown in Fig. 5. Although the general 




Figure 5: The comparison of pion form factors within rainbow-ladder truncation and with pion back- 
coupling included. 
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behavior of the form factor is already satisfying, we observe sizable deviations from one at zero 
momentum transfer. This violation of current conservation stems from the omission of the pion 
exchange as well as the seagull diagram in Fig. 4 and may indicate the size of corrections to be 
expected from these contributions. We are currently in the process of completing the calculation of 
the two remaining diagrams. In the future we expect to be able to extend the framework to other 
light- and heavy-meson channels as well as into the baryon sector. 

3. Tetraquarks 

3.1 Bethe-Salpeter equation for the tetraquark 

Bound states can be characterized as poles in the appropriate n-point Green function. The 
amplitude and the mass of a bound state are found by solving a Bethe-Salpeter equation [34]. The 
general structure of a multi-particle Bethe-Salpeter equation and especially the tetraquark BSE can 
be found in [35, 36]. They have the common structure: 

X ¥ = KG X ¥. (3.1) 

*F denotes the amplitude of the tetraquark bound-state, K is the four-body interaction kernel and Go 
represents a product of four fully-dressed quarks. The multiplication is defined in a functional sense 
and assumes implicit integration of all intermediate momenta and Dirac, color and flavor indices. 
Following the successful strategy in the baryon sector [37], we drop all 3PI and 4PI contributions. 
The resulting truncated kernel is a sum of the three remaining pair interactions [36]: 

K = Y^Kaa' and K aa , = K a + K a , - K a K a , . (3.2) 

aa 1 

The subscripts a and a' denote qq, qq and qq pairs, respectively. In principle Eq. (3.1) can now 
be solved with the same techniques used for the nucleon [38]. This is, however, a numerically 
demanding undertaking and therefore a reduction to a two-body system is employed. 

The construction of an appropriate 2-body system involves the pair-interacting scattering T- 
matrix 

Taa> = T a + T a , + T a T a , (3.3) 
which is closely related to the interaction kernels via Dyson's equation: 

7^=^(1 + 7^) and T a = K a (l + T a ), (3.4) 




Figure 6: Tetraquark BSE in the meson-meson/antidiquark-diquark picture. The hatched amplitudes involve 
diquark quantities; the remaining ones are of mesonic nature. Single (double, dashed) lines are dressed quark 
(diquark, meson) propagators. 
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where we drop the first two terms in Eq. (3.3). It is important to note that the two-body scattering 
matrix T a contains a pole of either a meson or diquark: 

T a = -r a D a r a , (3.5) 

where T is the meson or diquark amplitude and D the corresponding propagator. These pole con- 
tributions are assumed to be dominant. The key idea in the construction of a two-body equation is 
now the substitution of the interaction kernels with the T-matrices in a pole approximation. This 
ensures the BSE to be expressed in previously calculated degrees of freedom: mesons, diquarks 
and quarks. 

Putting everything together, the BSE in Eq. (3.1) reduces from a four-body equation, featuring 
quarks with gluons as exchange particles, to an effective two-body equation. The constituents are 
now mesons and diquarks interacting via quark exchange, see Fig. 6. We take into account only 
the mesons and diquarks with lowest mass, expecting the higher-mass state contributions to be 
subleading. It is interesting to note that our approach does not permit a pure diquark-antidiquark 
state. Combining both equations in Fig. 6 renders diquark-antidiquark contributions to appear 
internally only. Thus, one may view the resulting tetraquark bound state as a meson molecule with 
diquark-antidiquark admixture to its kernel. 

3.2 Results and Discussion 

Our result [39] for the mass of the up/down ++ tetraquark state as a function of the pseudo- 
scalar-meson mass is shown in the left panel of Fig. 7, together with a calculation that only includes 
the meson-molecule component of the tetraquark. The line represents a fit to the data including a 
constant, a square root and a linear term. Besides favoring a 7T7r-molecule picture over a sublead- 
ing diquark-antidiquark component, the main result of our investigation is the value for the u/d 
tetraquark at the physical point with nips = m n : 

m" T /d (0 ++ ) =403MeV. (3.6) 




ni ps [Gel/] m Quark iGeV] 



Figure 7: Mass of the ++ tetraquark as a function of the pseudoscalar-meson mass (left) and the quark 
mass [39]. The large red circle highlights a potential all-charm tetraquark. 
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This value is somewhat lower than the mass deduced from the Roy-equation [13, 14] but other 
bound states that can mix with the scalar tetraquark are not included in our work. The behavior of 
the tetraquark mass in dependence of the quark mass shows some interesting features: Except for 
the chiral mass region, the fit resembles the typical behavior of a Goldstone boson. This can be 
readily understood from the strongly dominating meson-meson component to the tetraquark bound 
state. Furthermore, the large decay constant and the low mass of the a is in accordance with our 
picture of a 7T7r-molecule dominated tetraquark. 

It is furthermore interesting to speculate about the existence of an all-charm tetraquark state. 
In Fig. 7 the largest pseudoscalar-meson mass corresponds to a quark mass in the charm region. 
We therefore read off the mass of an all-charm scalar tetraquark state to be at 

m^(0++)=5.3GeV. (3.7) 

This mass is considerably lower than the 6.2 GeV which were obtained in simpler model calcula- 
tions [40, 41]. It is also much lower than the r] c threshold. Potential decay channels into D mesons 
and pairs of light mesons necessarily involve internal gluon lines. This could result in rather small 
decay widths. Further results for tetraquarks with different quantum numbers are in progress. 
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